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We use novel STM techniques in concert to study the doping dependence of electronic structure 
in Bi2Sr2CaCu208+5. At all dopings, the low energy states are relatively homogenous except for 
dispersive density-of-states modulations whose properties are used to elucidate the momentum- 
space characteristics of quasi-particles. The superconductive coherence-peaks, ubiquitous in near- 
optimal tunneling spectra, are destroyed with strong underdoping. A new spectral type, likely 
characteristic of the zero temperature pseudogap regime, appears in these samples. Exclusively 
in regions exhibiting this new spectrum, we find quasi periodic modulations in dl/dV as well as 
in topography, with Q = (±27r/4.5ao, 0) and (0, ±27r/4.5ao) ± 15%. This is consistent with the 
existence of a local charge density modulation at these wave vectors. Surprisingly, this state coexists 
harmoniously with the low energy nodal quasi-particles. We discuss the relevance of these findings 
to the cuprate phase diagram and to the relationship between the pseudogap and superconductivity. 



High temperature superconductivity emerges in the 
cuprates when localized electrons in the Cu02 plane be- 
come itinerant due to hole-dopingi A schematic diagram 
of the dependence of their electronic/magnetic character- 
istics on temperature (T), and doped- hole density (p), is 
shown in 

Fig. 1 (adapted from Ref. 2). This is referred to col- 
loquially as the cuprate 'phase-diagram' although true 
electronic phases (antiferromagnetic insulator (AFI), su- 
perconductor (SC) and metal (M)) have been identi- 
fied in only three of its regions. Among the other 
poorly understood and sometimes overlapping regions 
are those characterized by a 'pseudogap' (PG); -^i^i'^ as 
a disordered magnetic 'spin glass' (SG) with incommen- 
surate spin fluctuations and as a non-Fermi liquid 
(NFL)ii^ Although no widely accepted theory exists for 
these phenomena, research into the underdoped cuprates 



{p < 0.16) is primarily motivated by two (apparently) 
distinct points of view. In the first, another electronic 
phase (ordered state) competes with the superconductiv- 
ity. This state is thought to drive the superconducting 
critical temperature Tc to zero as it strengthens with re- 
duced doping, and to disappear near maximum Tc at a 
hidden quantum critical point (X in Fig. 1). Within 
this picture, the pseudogap is primarily a property of 
an unknown electronic ordered state. The second pic- 
ture is one in which, because of low superfluid density at 
low doping in two dimensions, superconductivity is de- 
stroyed by order-parameter phase and/or amplitude fluc- 
tuations and, as superfluid density rises with increased 
hole-doping, the T^ rises. From this perspective, the 
pseudogap characterizes a region above Tc exhibiting 
gaps in the spin and charge spectrum due to the exis- 
tence of singlet pairs, even though the material is not 
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a superconductor. Here we report on doping-dependent 
STM studies of Bi2Sr2CaCu208+5(Bi-2212), designed to 
explore these two models. 

A serious complication for experimental exploration of 
these ideas is the fact that real-space (r-space) probes 
generally detect nanoscale spatial heterogeneity in the 
electronic/magnetic structure of underdoped cuprates. 
For example, when 0.03 < p < 0.14, muon spin rota- 
tion (^SR) studies indicate the existence of a disordered 
magnetic 'spin glass' in La2_KSr^Cu04 (La-214) and 
oxygen doped LaCu04-|-57 and in Bi2Sr2CaCu208+5(Bi- 
2212) 1^ Scanning tunneling microscopy (STM) based 
local-density-of-states (LDOS) imaging reveals nanoscale 
electronic structure variations in (Bi-2212)2ii2iii4i2ii^ and 
in Naa;Ca2-KCu02Cl2r-" Nuclear magnetic resonance 
(NMR) points to strong nanoscale carrier density dis- 
order with variations in local p of at least 25% of mean 
carrier density in both underdoped La-214i^*i& and un- 
derdoped Bi-2212J^ For a wide variety of underdoped 
cuprates, scaling analyses of penetration depth measure- 
ments reveal finite size effects consistent with nanoscale 
heterogeneity in the superfluid density^ These results 
provide abundant independent evidence that spin and 
charge degrees of freedom are heterogeneous at the 
nanoscale in many underdoped cuprates. It has not 
yet been possible to determine if this heterogeneity is 
a sample-specific and extrinsic effect due to crystal, 
dopant, or chemical disorder, or is an intrinsic effect of 
the cuprate electronic structure. Nor have its implica- 
tions for the phase diagram been considered widely. 

A complementary description of electronic structure 
to that in r-space is in momentum space (fc-space), 
accessible for cuprates via angle resolved photoemis- 
sion (ARPES), and optical techniques. ARPES reveals 
that, at optimal p in the superconducting phase, the 
Fermi-surface (FS) of hole-doped cuprates is gapped 
by an anisotropic energy gap A(fc) with four nodes, 
and, below Tc, quasiparticles exist everywhere along the 
normal-state Fermi surface li^iSSiSi At a fixed low tem- 
perature, those quasi-particle states with k — {ztn/aQ, 0) 
and (0, ±7r/ao) near the 1st Brillouin zone-face, degrade 
rapidly in coherence with reduced doping until they 
have become incoherent at p < O.lOii^iSSiSSiS^ This is 
a mysterious phenomenon, since it is closely correlated 
with superfluid density^i but, within conventional su- 
perconductivity theories, should not be related to it. 
By contrast, states on the 'Fermi-arc' (FA)2^ nearby 
the nodes retain their coherence down to the lowest 
dopings studied liSiSSiSi Transient grating optical spec- 
troscopy studies of non-equilibrium quasi-particles in un- 
derdoped YBa2Cu306.3^ also find hfetimes for antinodal 
excitations that are orders of magnitude shorter than 
those of the nodal quasi-particles. Thus, the electronic 
structure of underdoped cuprates also appears to be het- 
erogeneous in fc-space, in the sense that states proximate 
to the gap-nodes (nodal) have quite different character- 
istics and evolution with doping, than those near the 
zone- face (antinodal). No generally accepted explana- 



tion exists for either these fc-space phenomena or their 
relationship to the phase diagram. 

An electronic structure like this, which is heteroge- 
neous in both r-space and fc-space, presents special chal- 
lenges to conventional experimental techniques. Never- 
theless for the cuprates, it is critical to understand the 
relationship between these two electronic structure repre- 
sentations and how it evolves with doping. Here, by using 
several innovative STM techniques in concert, we have 
carried out the first comprehensive study of the doping 
dependence of Bi-2212 superconducting electronic struc- 
ture in r-space and determined its relationship to the 
electronic structure in fc-space. Our first technique is lo- 
cal density of states {LDOS) mapping which consists of 
measurement of the tip-sample differential tunneling con- 
ductance g{r,V) = ^|r,v at each spatial location and at 
each sample bias voltage V. Since LDOS{f, E = eV) cx 
g{r,E = eV), an energy-resolved r-space image of the 
electronic structure is attained. Second, from g{r, E) a 
dataset, the magnitude of the energy-gap in the density of 
states A(r) can be determined in a process called a gap- 
map. Here A(r) = and is the energy of the 
first maximum in LDOS above(below) the Fermi level 
(neglecting impurity states)^ A final technique, recently 
introduced to cuprate studiesjSS is Fourier transform 
scanning tunneling spectroscopy (FT-STS). Here, the q- 
vectors of spatial modulations in ^(r, E) are determined 
from the locations of peaks in g{q^ E), the Fourier trans- 
form magnitude of g{f,E). This technique has proven 
valuable by virtue of its exceptional capability22i2ii2S to 
relate the nanoscale r-space electronic structure to that 
in fc-space. 

To clarify these techniques, we show examples of each 
measurement type in Fig. 2. All are derived from a single 
g{f, E) data set measured with 1.8 Aresolution in 256^ 
pixels on a 50 nm -square field of view (FOV). Figure 
2A is the gapmap A(r) from this g{f,E) using a color 
scale that spans 20meV < A(r) < 70 meV . The super- 
posed inset shows the resolution of a topographic image 
taken precisely where this g{f,E) was measured. An ex- 
ample of unprocessed g{r, E) typical of strongly under- 
doped samples, is shown in Fig.2B. These spectra were 
measured along the red line in Fig 2A. From these data, 
one can see dramatic changes during which the coher- 
ence peaks in the spectrum disappear and are replaced 
by to a distinctly different type of spectrum. Such evo- 
lutions are ubiquitous in strongly underdoped Bi-2212 
samples. Fig.2C is the measured gifjE = —12 meV) in 
the FOV of 2A. It exhibits a complex pattern of very 
weak LDOS modulations which are relatively homoge- 
nous in the sense that RMS-deviations from spatially av- 
eraged value of g{f, E) are only 10%. Fig. 2D is the 
g{q,E = — 12 meV) calculated from Fig.2C and it re- 
veals that the LDOS'-modulations in Fig.2C were made 
up of only a small set of well-defined g- vectors (^2*^ 

We apply these techniques in concert to study the 
doping-dependence of electronic structure in a series of 
Bi-2212 samples. They are all single crystals grown by 
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the floating zone method. Doping is controlled by oxy- 
gen depletion so that no other elemental impurities are 
introduced. Each is cleaved in cryogenic ultra-high vac- 
uum before immediate insertion in to the STM head. If 
its BiO surface is flat and free of nanoscale debris, each 
sample is usually studied for several months, typically 
in a 50nm square field of view (FOV). More than 10^ 
spectra were acquired for the studies reported here. 

In Fig. 3 we show 50nm-square gapmaps measured on 
samples with five different dopings. Identical color scales 
representing 20 meV< A(r) < 70 meV are used for all 
images. The local hole concentration is impossible to 
determine directly, but we estimate that the bulk dop- 
ings were approximately 3A(0.19±0.01), 3B(0.18±0.01), 
3C(0.15 ± 0.01), 3D(0.13 ± 0.01), 3E(0.11 ± 0.01). Near 
optimal doping (Fig 3A,B) the gapmaps are heteroge- 
neous but nonetheless the vast majority of tunneling 
spectra are manifestly those of a superconductor (see be- 
low). However, at the lowest dopings and for gap values 
exceeding approximately 65 meV, there are very many 
spectra where A actually becomes ill defined because co- 
herence peaks do not exist at the gap edge (see for exam- 
ple Fig.3F, spectrum 6). We represent these spectra by 
black in the gapmap, since they are almost identical to 
each other and appear to be the limiting class of spectra 
at our lowest dopings. 

The spatially averaged value of A(r) for each crys- 
tal, A, and its full width at half maximum, cr, are: 
3A(A = 33 ± 1 meV, cr=7 meV), 3B(A = 36 ± 1 meV, 
cr=8 meV), 3C(A = 43 ± 1 meV, cr=9 meV), 3D( 48 ± 1 
meV, cr=10 meV), and 3B(A > 62 meV but with a ill de- 
fined). As doping is reduced, A grows steadily consistent 
with other spectroscopic techniques, such as ARPES and 
break-junction tunneling ji2i22i2^ which average over the 
heterogeneous nanoscale phenomena. This observation is 
very important because it demonstrates that our Bi-2212 
surfaces evolve with doping in an electronically equiva- 
lent fashion to those studied by the other techniques, and 
that we are probing the low temperature state of the un- 
derdoped pseudogap regime. 

In Fig.3F we show a series of the 'gap-averaged' spec- 
tra. Each is the average spectrum of all regions exhibiting 
a given local gap value (from the single 50nm FOV of 
Fig.2C). They are color-coded so that each gap-averaged 
spectrum can be associated with regions of the same 
colour in all gapmaps (Fig.'s 3A-E). The spectra are la- 
beled from 1-6, with numbers 1 through 4 providing clear 
examples of what we refer to as coherence peaks at the 
gap edge (indicated by the arrows). These gap-averaged 
spectra are consistent with data reported previously by 
from gapmap studies by Matsuda et ali^ Here, from our 
doping dependence study, we can report that this set of 
gap-averaged spectra is almost identical for all dopings. 
The dramatic changes with doping seen in A(r) (Fig. 3) 
occur because the probability of observing a given type of 
spectrum in Fig.3F evolves rapidly with doping. For ex- 
ample, the gap-averaged spectrum labeled as 1 in Fig.3r 
has a 30% probability of occurring in gapmap 3 A, 25% 



in 3B, 5% in 3C, less than 1% in 3D, and 0% in 3E. The 
spectrum labeled 6 has a 0% probability of occurring in 
3A, 0.1% in 3B, 1% in 3C, 8% in 3D and > 55% prob- 
ability in 3E. The evolution of these gapmaps (Fig. 3) 
with falling doping, from domination by heterogeneous 
but predominantly superconducting spectral characteris- 
tics (Fig.3A,B) to domination by spectra of a very dif- 
ferent type (Fig.3E) is striking. 

Despite the intense changes with doping in the 
gapmaps (whose information content is, by definition, de- 
rived from the coherence peaks a.t E — A(r)), the LDOS 
at energies below about 0.5A remains relatively homoge- 
nous for all dopings studied. Figure 3F reveals this low- 
energy LDOS homogeneity because, independent of gap 
value, the g{E) below w 25 meV are almost the same 
everywhere and for all spectra. These low energy LDOS 
do, however, exhibit numerous weak, incommensurate, 
energy-dispersive, spatial LDOS'-modulations with long 
correlation lengths (for example Fig.2C). We focus on 
the doping dependence of these low energy g(f, E) data 
by applying the FT-STS technique. Figure 

4A-C shows measured g{q^ E) for the three g{r^ E) 
datasets used to generate Fig.3A, 3D, and 3E. Each sub- 
panel is the measured g{q, E) at the labeled energy, with 
the reciprocal space locations of the Bi (or Cu) atoms 
q — (±27r/ao,0) and (0, ±27r/ao), appearing as the four 
dark spots at the corners of a square. It is obvious that 
multiple sets of dispersive LZ30S'-modulations exist at 
all three dopings, but each exhibits different trajectories 
as a function of E for different p. 

Analysis of these low energy LDOS modulations 
requires a model for their relationship to states in 
/c-space. We apply the "octet model" of quasipar- 
ticle interference^SiSi which is predicated on a Bi- 
2212 superconducting band-structure exhibiting four sets 
of 'banana'-shaped contours of constant quasiparticle- 
energy surrounding the gap nodesiiS. Because of the 
quasiparticle density of states at E is 

n{E) oc / ^-^dk (1) 

JE{k)=uj VkE{k) 

while each 'banana' exhibits its largest \'S/ kE{k)\~^ near 
its two ends, the primary contributions to come from 
the octet of momentum-space regions at the ends of each 
'banana' kj] j=l,2,..8 (Fig.5A). Mixing of quasiparticle 
states in the octet by disorder scattering produces quasi- 
particle interference patterns which are manifest as spa- 
tial LDOS'-modulations. The intensity of such scattering 
induced modulations is primarily governed by joint den- 
sity of states (among other factors). The wavevectors of 
the most intense LDOS-modulations are then determined 
by all possible pairs of points in the octet kj. Sixteen 
distinct -\-q and —q pairs should be detectable at each 
non-zero energy by FT-STS. From them, the energy de- 
pendence of the octet locations kj (E) can be determined 
and associated with a 'locus of scattering' ks{E). Com- 
prehensive internally-consistent agreement between Bi- 
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2212 STM data and this model is achieved near optimal 
doping. '^^ Until this work, nothing was known about its 
utility for strongly underdoped cuprates. 

Theoretical analyses beyond the simple octet 
model34i35,36,37,38,39,40,4i,42,43 capture many elements of 

our previously reported g{q,E) data, but no resolution 
of the exact source, strength, or type of scattering 
has yet been achieved. Nevertheless, the existence of 
numerous sets of long-correlation length, dispersive, 
LDOS modulations, all of which are self-consistent 
with a single A(fc) for both filled and empty states, is 
indicative of good Bogoliubov-like quasi-particles. Since 
the LDOS-modulations can be associated consistently 
with a 'locus of scattering' ks{E) via the octet model, 
we analyze our observations within this model using the 
g*- vector designations shown in Fig.5A. 

Figure 5B shows the measured length of qi, 55 and qy 
as a function of energy for the three datasets in Fig. 4. 
Figure 5C shows the locus of scattering calculated for 
these three using: 

qi = (2fc^, 0); 95 = (0, 2ky);qi = {k^ - ky, ky - k^) (2) 
fc, = {±k^{E), ±ky{E))- ks = i±ky{E),±k4E)) (3) 

The ks{E) determined by this technique differs only 
slightly between dopings. even though the actual g(f, E) 
for different dopings are quite different at any given en- 
ergy. These three ks{E) are each the same for filled and 
empty sates and each consistent with the same A(fc) at 
that doping. Thus Bogoliubov-like quasi-particles appear 
to exist at these momentum space locations at all dop- 
ings. This is consistent with the small motion of the FS 
in this doping range detected by ARPES. 

These observations certainly do not exhaust the 
changes observed in g{f,E) with falling doping. A very 
strong effect is the evolution, with doping, of the q- 
space location of strongest LZJOS'-modulation at any en- 
ergy. This modulation is always associated with qj and 
vits location evolves from qi = ^ at p= 0.19 ± 0.01, 
to 9i = 51^ at p= 0.14 ± 0.o7, to <fi = 4^ at 
p= 0.10 ± 0.01. Another effect is a decrease in relative 
intensity of dispersive LDOS'-modulations (72, 93, 96, 97 
relative to those of qi, 95, with falling with p. These 
effects will be reported in detail elsewhere. 

The doping dependence of states with k — (±7r/ao, 0) 
and (0, ±7r/ao) in the 'flat band' region near the zone- 
facsiS (green shaded areas in Fig. 5 A) is extremely dif- 
ferent. These states can also be identified by FT- 
STS analysis of g(r, E) data. By definition, the co- 
herence peaks in g(f,E) occur aX E = A(r). In 
all samples, they exhibit intense particle-hole sym- 
metric LDOS'-modulations, with wavevectors G = 
(±27r/ao,0) and (0, ±27r/ao)i^ These coherence peak 
LDOS'-modulations at i? = A possibly occur due 
to Umpklapp scattering between k = (±7r/ao,0) and 
(0,±7r/ao)^. Therefore, the coherence peaks in tunnel- 
ing are identified empirically with the zone-face states 



at A: = (±7r/ao,0) and (0,±7r/ao). This identification 
is also consistent with theory. The coherence peaked 
tunneling spectra (e.g. Fig.2F: spectra 1-4) are theo- 
retically viewed as due to superconducting pairing on 
the whole FS because such spectra are consistent with a 
Aj.2_j,2 everywhere on the ARPES-determined FS near 
optimal doping.'' Wc therefore consider any spatial re- 
gions that show clear coherence peaks with q — G LDOS- 
modulations, to be occupied by a canonical d-wave su- 
perconductor (dSC). 

Near optimal doping, more than 98% of any FOV ex- 
hibits this type of coherence peaked dSC spectrum. As 
the range of local values of A(r) increases with decreas- 
ing doping, the intensity of the q = G coherence peak 
LDOS'-modulations becomes steadily weaker until, wher- 
ever A(r) > 65 meV, they disappear altogether. This 
process can be seen clearly in the gap-averaged spec- 
tra of Fig.SF where the average height of the coher- 
ence peaks declines steadily with increasing A. It is 
found equally true for all dopings. Wherever the co- 
herence peaks and their q = G i 13 05*- modulations are 
absent, a well-defined new type of spectrum is always ob- 
served. Figure 6A shows a high resolution gapmap from 
a strongly underdoped sample. Examples of this new 
type of spectrum, along with those of coherence peaked 
dSC spectra, are shown in Fig.6B. The coherence peaked 
spectra (red) are manifestly distinct from the novel spec- 
tra (black) which have a V-shaped gap reaching up to 
-300 meV and +75 meV. For reasons to be discussed be- 
low, we refer to the new spectrum as the zero temperature 
pseudogap (ZTPG) spectrum. 

The replacement of coherence peaked dSC spectra by 
ZTPG spectra first begins to have strong impact on av- 
eraged properties of g{r, E) and g(g, E) below about 
p=0.14 where the fractional area covered by ZTPG 
spectra first exceeds 10% of the FOV. In terms of 
the spectral shape no further evolution in the form of 
the ZTPG spectrum is detected at lower dopings. In- 
stead, a steadily increasing fractional coverage of the 
surface by these ZTPG spectra is observed. Our pre- 
vious studiesiSii^ were carried out at dopings p > 0.14 
and, when ZTPG spectra have been detected at such 
higher dopings it is in a tiny fraction of the FOV. 
Very significantly, spectra similar to the ZTPG spectrum 
are detected inside cores of Bi-2212 quantized vortices 
where superconductivity is destroyed .^4^^ Furthermore, 
a very similar spectrum is observed in another very un- 
derdoped cuprate Na^Ca2_a;Cu02Cl2r^*^ even in the 
non-superconducting phase. It therefore seems reason- 
able that the characteristic zero-temperature spectrum 
in the pseudogap phase is of this type. This is why we 
tentatively assign this spectrum the ZTPG designation. 

As discussed above, the g{f,E) for E < A/2 in even 
the most underdoped samples (Fig.SE, Fig.4C) exhibit 
relatively homogenous electronic structure with good 
quasi-particles dispersing on the Fermi-arc. However, for 
E > A/2 in these same samples, our previous analyses 
techniques fail, probably because very different phenom- 
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ena are occurring in different nanoscale regions of each 
FOV. To explore the implications of the ZTPG spec- 
trum for strongly underdoped samples, new analysis tech- 
niques are therefore required. Here we introduce a mask- 
ing process which has proven highly effective. From a 
given strongly underdoped data set, the g{r, E) in all re- 
gions where -EA > 65 meV are excised and used to form a 
new masked data set g{r, -E')1a>65- The remainder forms 
a second new dataset 5((r, i?)|A<65- The > 65 meV 
cutoff was chosen because, on the average, it represents 
where the coherence peaks with associated q = G mod- 
ulations have all disappeared and are replaced by the 
ZTPG spectra. An example of this type of mask for the 
gapmap in Fig.6A is show in Fig.5C. It is important to 
note a serious drawback of the masking process. The 
q resolution is of masked data is considerably worse in 
than those shown in Fig.5B because the largest contigu- 
ous nano region in the mask is about 20% of the full 
FOV. As a result, the precise modulation period and dis- 
persions of any effects detected by masking cannot can 
be determined with nearly the same degree of accuracy 
as the low-energy quasi-particle interference signal (Fig. 5 
and Ref. 28). 

FT-STS analysis of such (^(r, £;)|a<65, 5(^,^)|a>65) 
pairs shows that they exhibit dramatically different phe- 
nomena. In the -E)|a<65, the dispersive trajectory 
of qi is seen up to « 36 meV and no further LDOS- 
modulations can be detected at any higher energy (red 
symbols in Fig.6D). In the g(r, i?)|A>65 data, the identi- 
cal dispersive qi signal is seen below E 36 meV. However, 
a new LUOS'-modulation appears in the g(r, -E)|a>65 be- 
tween E> 65 meV and our maximum energy E=150 meV 
(black symbols in Fig.6D). We designate its wavevector 

r- 

To explore the real-space structure of this new high- 
energy LZJOS'-modulation, we define a map rg|*^(r) — 
^£=655(^' ^)l<5>65 which sums over this energy range. 
This map is shown in Fig.7A and, although it has quite 
a disordered mask, careful examination reveals checker- 
board modulations within each nano region. Impor- 
tantly, Fourier transform analysis of this rg5°(r) shown 
in Fig. 7B reveals a well-defined wavevector set q* — 
(±27r/4.5ao,0) and (0, ±27r/4.5ao) ± 15% for these new 
high-energy modulations as indicated by the arrow on 
black data points in inset to Fig. 7B. The identical anal- 
ysis the complementary map Sg|°g(r, i?)|5<65 is feature- 
less near q* (red data in inset of Fig. 7B). Thus, an 
LZ305'-modulation with very low (or zero) dispersion, 
exists above ±65 meV exclusively in regions character- 
ized by ZTPG spectra of strongly underdoped Bi-2212 
samples. 

Constant-current topography represents, albeit loga- 
rithmically, the contour of constant integrated density 
of states up to the sample-bias energy. It does not 
suffer from the systematic problems due to effects of 
the constant-current setup condition renormalizationiS 
which plague g{r,E). It is therefore a more conclusive 
technique for detection of net charge density modulations 



by STM. To search for topographic modulations, we ap- 
ply the identical mask (Fig.7C) to the topographic image 
which was acquired simultaneously with the gapmap in 
6 A. The magnitude of the Fourier transform along the 
g||(27r,0) for the masked topographic image shows that, 
in the ZTPG regions, the topography is modulated with 
qtopo = (±27r/5ao,0) and (0,±27r/5ao) ± 25%(indicated 
by the arrows in Fig.7D). Fourier transforms of the 
complementary part of the topography (from A <65 
meV regions not exhibiting ZBTG spectra) shows no 
such modulations at any wavelengths near this qtopo (red 
in Fig.7C). No topographic modulations near qtopo are 
detected anywhere in samples at higher doping. An 
additional peak near qrecon — (27r/2ao, ±27r/2ao) and 
(— 27r/2ao, — 27r/2ao) in the Fourier transform of the to- 
pograph comes from a reconstruction along the super- 
modulation maximum, is observed at all dopings and, 
since no signature is observed in LDOS at qrecon for any 
energy or doping, is regarded as irrelevant. 

Some important new facts about the strongly under- 
doped regime of Bi-2212 electronic structure emerge from 
these results. Our first finding is related to the Fermi- 
arc quasi-particle states. As shown in Fig. 5C, FT-STS 
indicates that quasi-particle interference occurs between 
Bogoliubov-like states in approximately the same region 
of /c-space for all dopings. These Fermi-arc quasiparti- 
cles remain spatially homogenous (except for relatively 
weak LDOS-modulations) even in the most underdoped 
samples. They are Bogoliubov-like in the sense that 
they exhibit particle-hole symmetry at each location in 
fc-space, and, at each doping, are all consistent with the 
same A(fc). Therefore, one can reasonably postulate that 
Fermi-arc states are gapped by superconducting interac- 
tions at all dopings studied. If so, they, and the asso- 
ciated 'nodal' superconducting state^^iiSiiSiffi are amaz- 
ingly robust against the heterogeneous electronic phe- 
nomena which so dominate Bi-2212 at other energies. 

Our second finding is the very different fate of states 
in the flat-band regions near k « (±7r/ao,0); (0, ±7r/ao). 
The appearance of ZTPG spectra in strongly underdoped 
samples coincides exactly with destruction of antinodal 
superconducting coherence peaks. Exclusively in these 
ZTPG regions (black in the gapmap), we find the three 
new modulation phenomena: (1) topographic modula- 
tions with qtopo = (±27r/5ao,0) and (0, ±27r/5ao) ±25%, 
(2) a peak in LDOS centered around q* = (±27r/4.5ao, 0) 
and (0,±27r/4.5ao) ± 15% for E >65 meV to at least 
E=150 meV and, (3) the dispersive qi quasi-particle 
branch exhibits its maximum modulation intensity when 
it passes through qi — qtopo = q* (Fig 7D). This last sit- 
uation has been predicte d'^^i'^^i^^i'^^i^^ as a consequence of 
the Fermi surface geometry and quasi-particle dispersion 
in the presence of potential scattering from charge-order 
with fixed Q; in this case Q = qtopo — q*- Taken together, 
these observations all point to the appearance of an un- 
usual charge ordered state with q* = (±27r/4.5ao, 0) and 
(0, ±27r/4.5ao) ± 15%, occurring only in the regions char- 
acterized by the ZTPG spectrum and only in strongly 
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underdoped Bi-2212. Note that due to the strong disor- 
der and the hmited size of the mask domain we can not 
distinguish between a charge density modulation caused 
by a condensed charge order and that caused by impu- 
rity through an enhanced charge susceptibihty in a state 
without condensed charge order. In addition we do not 
imply that there is true charge long range order since 
from the Imry-Ma argument this is always absent in a 
disordered system like Bi-2212. 

Charge order has been observed in other under- 
doped cuprates, including Nd doped La-214 with inelas- 
tic neutron scattering^ and more recently by STM in 
Naa;Ca2-a:Cu02Cl2i^ It has also been proposed, based 
on reported nondispersive (between and 20 meV), 'line 
object', LDOS modulations with q = (27r/4ao, 0)^^, that 
static stripes exist in optimally doped Bi-2212 below Tc. 
However, none of these phenomena have been detected 
in several independent higher resolution studiesi^fi'^i'^ 

Other very suggestive findings have also been made 
by STM in Bi-2212. Field induced sub-gap LDOS- 
modulations, with qvortex = (±27r/4.3ao, 0) and 
(0, ±27r/4.3ao)±15% were discovered surrounding vortex 
cores (where superconductivity is destroyed) near opti- 
mal doping. This observation provided the first STM 
evidence for some type of incipient charge-order com- 
peting with superconductivity in cupratesiSS Pioneering 
STM experiments to map the low energy LDOS above Tc 
when superconductivity is also destroyed, have detected 
sub-gap LDOS-modulations with qpq = (±27r/4.6ao, 0) 
and (0, ±27r/4.6ao)^ at near-optimal doping. Although 
neither of these low energy phenomena are fully under- 
stood, they do, along with low energy phenomena in 
ZTPG regions reported here, form a triad of apparently 
consistent observations. Destruction of superconductiv- 
ity, whether by high magnetic fields, by exceeding Tc, or 
by strong underdoping, results in very similar effects on 
low energy LDOS modulations. It remains to be deter- 
mined how the vortex core and pseudogap observations 
relate directly to the charge order. 

The identity of the electronic phase represented by the 
Qtopo = Q* charge order is difficult to discern. In the 
absence of disorder, plaquette orbital-order phases such 
as staggered flux phase (SFP)^^*^ and D-density-wave 
(DDW),^^ or intra-plaquette orbital phasesjS are not 
expected to exhibit topographic or LDOS'-modulations. 
However, in theory, LZJOS'-modulations can be pro- 
duced by vortex and disorder scattering in the SFP 
and DDW phases but not near q = (±27r/4ao,0) and 
(0, ±27r/4ao).'i2iiiS It remains theoretically unexplored 
whether disordered orbital phases could result in the 
complete set of new phenomena (see below) we report 
here. Charge-ordered phases including stripeSfSMiiS^ dis- 
order pinned electronic liquid crystal^ strong-coupling 
spin- and charge-density waves^i and, recently, hole-pair 
crystals^ have been proposed to exist in underdoped 
cuprates. Each of these phenomena would yield both 
LDOS- and topographic-modulations. But again, it re- 
mains theoretically unexplored whether these theories 



can account for our observations that in the ZTPG re- 
gions (i) a characteristic new tunneling spectrum exists, 
(ii) the q* modulations appear only above a relative high 
energy (w 65 meV) , (iii) they exhibit an incommensurate 
wavevectorq = (±27r/4.5ao, 0) and (0, ±27r/4.5ao)±15%, 
(iv) they exhibit the same spatial phase for positive and 
negative biases - so that the filled-state density maxima 
coincide with the empty-state-density maxima, and (iv) 
the q* modulations are replaced by the dispersive quasi- 
particle interference signals at sub-gap energies. A fur- 
ther point is that the predicted strong breaking of the 90° 
rotation symmetry in the stripe scenario is not observed 
in any of the STM studies, but again, it may be possible 
that this is due to the presence of strong disorder. For 
all these reasons the precise identity of the charge-order 
state in strongly underdoping Bi-2212 remains elusive. 

The data reported here also motivate a new conjec- 
ture on the evolution of electronic structure with reduced 
doping in Bi-2212. In r-space, we identify two extreme 
types of LDOS spectra (Fig. 6B) . The first exhibits clear 
coherence peaks at the gap edge and dominates near- 
optimal samples. The second type (ZTPG-spectrum) ex- 
hibits a V-shaped gap over a much wider energy range 
and dominates in strongly underdoped samples. We as- 
sociate the former with a pure d-wave superconducting 
state and conjecture that the latter reflects a zero tem- 
perature charge ordered state existing at sufficiently low 
dopings in the pseudogap regime. If the whole Bi-2212 
sample were homogeneous and consisted of only one of 
the above phases, then, in fc-space, quasi- particle peaks 
would exist all along the Fermi surface in the pure d-SC 
phase, but only on a finite arc around the gap nodes in 
the charge ordered phase with the zone face states be- 
ing incoherent due to localization. This may indeed be 
the case in Na:rCa2-:ECu02Cl2.2L^ In Bi-2212 he real- 
ity is more complicated. In optimally doped samples, 
more than 98% of the surface area exhibit dSC spectra 
(Fig. 3A,B). As doping falls the ZTPG regions appear 
and grow in significance until as doping approaches p 0.1, 
almost 60% of the area exhibits the ZTPG characteris- 
tics and the associated charge-order (Fig. 3E). From this 
trend, it is reasonable to expect that, at even lower dop- 
ing in the zero temperature pseudogap phase, 100% of 
the sample would exhibit the ZTPG characteristics and 
be charge-ordered. 

In a spatially disordered situation, the probability of 
occurrence of the two types of phenomena evolves con- 
tinuously with doping in a fashion related to the evo- 
lution of the gapmaps in Fig. 3. Therefore, properties 
which average over nanoscale phenomena would appear 
to evolve smoothly between the two extremes. Experi- 
mental results which are relevant to this proposal include, 
for example, the doping dependence of the ARPES (tt, 0) 
pgg^]j^2ji9j20j24 ^jjg specific heat jump at the superconduct- 
ing phase transition,™ and the c-axis conductivityj^ Due 
to the heterogeneous mixture of the dSC and the ZTPG 
regions, it is difficult for a spatially averaged experiment 
like ARPES to discern the properties of the ZTPG re- 
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gion. However, the well-known tendency of the coherent 
quasiparticle peaks near the zone face to be suppressed 
by underdopingfi2i22iS4 is consistent with our conjecture. 
The specific heat jump at the superconducting transition 
is a characteristic specifically of a dSC phase and not of 
the charge-ordered ZTPG phase. Hence, the declining 
specific heat jump as a function of underdoping^'^^ is 
also consistent with our hypothesis. Finally, due to tun- 
neling matrix element effects, c-axis tunneling senses the 
zone-face quasiparticles instead of the nodal ones. Hence 
the decrease of the c-axis conductivity with underdoping^ 
also seems consistent with our suggestion. Further inter- 
comparison between the results of these experimental 
techniques will be required to explore these proposals. 

One may also wonder if the charge-ordered ZTPG 
phase could actually be the source of the high tem- 
perature pseudogap phenomenal^ Since the charge-order 
would melt as a function increasing temperature, the an- 
swer to this question it is not clear. If the charge order 
melting temperature were actually the pseudogap tem- 
perature T*, then of course, the entire pseudogap phe- 
nomena would be caused by the ZTPG phase. On the 
other hand if the charge melting temperature were lower 
than the pseudogap temperature T*, then there would 
exist yet another mechanism that gaps the spin excita- 
tions at the pseudogap crossover. 

Independent of these conjectures on doping- 
dependence of spatially-averaged observables, the 
data reported here represent two significant advances. 
First, a charge ordered state competing with super- 
conductivity exists in strongly underdoped Bi-2212 
and second, the charge-ordered regions and the d-SC 
regions share the same low energy nodal quasiparticles. 
New microscopic theories for the electronic structure of 
underdoped cuprates are therefore required to explain, 
not only the unusual characteristics of the charge order 
phenomenon and how it competes with superconductiv- 
ity for spectral density at the zone-face, but also how it 
meshes so smoothly with the nodal superconductivity at 
low energies. 

Figure 1 

A schematic phase diagram of the cuprates. Note that 
the region of our studies is between 0.1 < p < 0.2 and 
near the zero temperature axis of this diagram where the 
magnetically disordered spin glass region coexists with 
the superconductivity. 
Figure 2 

A The gapmap A(r) calculated from a g{f,E) data set 
which was measured in a 55 nm field of view (FOV) with 
256^ pixels; inset topography in this FOV. B The un- 
processed linecut (;(r, E) connection a region with sharp 
coherence peaks with a region designated ZTPG (see 
text), along the red line in Fig. 2 A. C The measured 
g(r, E = — 12meV) in the identical FOV as 2A. D calcu- 
lated from C (the reciprocal space location of the Bi or 
Cu atoms are labeled (27r,0). These five types of mea- 
surements are used in concert to study the doping depen- 
dence of electronic structure in Bi-2212. 



Figure 3 

A-E Measured A(r) for five different hole-doping levels. 
F. The average spectrum associated with each gap value 
in a given FOV. They were extracted from the g{f, E) 
that yielded Fig. 3C but the equivalent analysis for 
g{f, E) at all dopings yields results which are indistin- 
guishable. The coherence peaks can be detected in #'s 
1-4. 

Figure 4 

A-C Examples of measured g{q, E) for a variety of ener- 
gies E as shown at three doping levels. 
Figure 5 

A. A schematic representation of the 1st Brillouin zone 
and Fermi surface location of Bi-2212. The flat-band re- 
gions near the zone face are shaded in blue. The eight 
locations which determine the scattering within the octet 
model are show as red circles and the scattering vectors 
which connect these locations are show as arrows labeled 
by the designation of each scattering vector. B. Mea- 
sured dispersions of the LlJOS'-modulations (fi, and 
qj for the 3 dopings whose unprocessed data is shown 
in Fig. 4. We chooses this set of three g-vectors because 
they exhibit the maximum intensity of any set sufficient 
to independently determine the locus of scattering ks [E) 
for all dopings. C. Calculated loci of scattering ks{E) for 
all 3 dopings. The blue line is a fit to the 89KOD data 
Figure 6 

A A high-resolution 46nm square gapmap from a 
strongly underdoped sample. B. Examples of represen- 
tative spectra from (1) nanoscale regions exhibiting co- 
herence peaked spectra with q — Q LDOS'-modulations 
aX E = A(f^ (red) and (2) from regions exhibiting ZTPG 
spectra (black). The locations where these spectra occur 
are shows as small dots in A. C. The mask identifying re- 
gions with A <65 meV from the ZTPG regions calculated 
from gapmap in A. D. Dispersion of qi{E) in regions 
with coherence peaked spectra A < 65 meV is shown in 
red. There are no modulations at any higher energies 
within our range. Dispersion of qi{E) in regions with 
ZTPG spectra for E< 36 meV (black symbols). The red 
squares in Fig.SB represent a combination of these two 
(indistinguishable) dispersions in this range. For E> 65 
meV, the wavevector of the new modulations in ZTPG 
regions are shown in black. To within our uncertainty 
they do not disperse and exhibit = (±27r/4.5ao, 0) 
and (0, ±27r/4.5ao) ± 15%. 
Figure 7 

AThe image of g{f, E) masked by Fig.6C and then 
summed from 65 meV to 150 meV ■T)^{r) = 
^65°3(^i -^)|i5>65- The light grey regions are outside the 
mask. This FOV is 46nm square and a careful ex- 
amination reveals a 'checkerboard' modulation occurring 
within the regions inside the mask. The equivalent im- 
age for r|5°(r) = T,lf^g{f, E)\s<:q5 is featureless except 
for Bi atom locations, and the supermodulation which is 
at 45° to the modulation in A. B. Fourier transform A. 
The Bi atom locations are circled in orange. A square 
of diffuse maxima are observed surrounding the central 
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point. The inset shows the plot of the Fourier transform 
amphtude along the line shown. It reveals a maximum 
at q* ^ (±27r/4.5ao,0) and (0, ±27r/4.5ao) ± 15% as in- 
dicated by the arrow. C. The magnitude of the Fourier 
transform of the masked topographic image, taken at 150 
mV and 150 pA, along the gtopo ||(27r, 0) direction is mod- 
ulated with qtopo = (±27r/5ao,0) and (0, ±27r/5ao) ±25% 
(black squares the black line is a guide to the eye). 
Fourier transforms of the complementary part of the to- 
pography (from < 65 meV regions not exhibiting ZBTG 
spectra) shows no such modulations at any wavelengths 
near this qtopo (red triangles the red line is a guide to the 
eye). The difference between these two Fourier transform 
intensities is shown in blue and shows the degree to which 
the ZTPG region shows topographic modulations unde- 
tectable elsewhere. The relative weakness of the topo- 
graphic modulations at qtopo is because of the logarithmic 



sensitivity of constant-current topography to net charge 
modulation. D. A plot of the amplitude of the qi LDOS 
modulation as a function of \qi \ for the sample data as in 
Fig.SC, 4E, 5D. The maximum intensity of the modula- 
tions in the ZTPG regions occurs at \qi \ = 27r/4.8aolO%. 
No special scattering is observed of the quasiparticles in 
the dSC regions. 
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